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ABSTRACT 

In an earlier investigation, we proposed population boundaries for both inspiralling and mass- 
transferring double white dwarf (DWD) systems in the distance independent "absolute" amplitude- 
frequency domain of the proposed space-based gravitational-wave (GW) detector, LISA. The de- 
generate zero temperature mass-radius (M-R) relationship of individual white dwarf stars that we 
assumed, in combination with the constraints imposed by Roche geometries, permits us to identify 
five key population boundaries for D WD systems in various ph ases of evolution. Here we use the non- 



zero entropy donor M-R relations of iDelove fc BildstenI (j2003l ) to modify these boundaries for both 



DWD and neutron star-white dwarf (NSWD) binary systems. We find that the mass-transferring 
systems occupy a larger fraction of space in "absolute" amplitude-frequency domain compared to 
the simpler T = donor mod el. We also di s cuss h ow these boundaries are modified with the new 



evolutionary phases found bv iDelove et al.l l|2007f l. In the initial contact phase, we find that the 



contact boundaries, which are the result of end of inspiral evolution, would have some width, as 
opposed to an abrupt cut-off described in our earlier T = model. This will cause an overlap be- 
tween a DWDs & NSWDs evolutionary trajectories, making them indistinguishable with only LISA 
observations within this region. In the cooling phase of the donor, which follows after the adiabatic 
donor evolution, the radius contracts, mass-transfer rate drops and slows down the orbital period 
evolution. Depending upon the entropy of the donor, these systems may then lie inside the fully 
degenerate T = boundaries, but LISA may be unable to detect these systems as they might be 
below the sensitivity limit or within the unresolved DWD background noise. We assess the limits 
and applicability of our theoretical population boundaries with respect to observations and find 
that a measurement of / by LISA at high frequencies (Log [/] > 2) would likely distinguish between 
DWD/NSWD binary. For low frequency sources, GW observations alone would unlikely tell us about 
the binary components, without the help of electromagnetic observations. 

Subject headings: accretion, accretion disks - binaries: close — gravitational waves — stars: white 
dwarfs, neutron stars 



Introduction 



Th e proposed space-based Gravitational-wave (GW) detector, LIs M {Laser Interferometer Space An 



Bender 



19981) . is sensitive to GWs in the 10" 



tenna) (jFaller fc Bendeilll984 lEvans. Iben. fc Smarilll987 : 

Hz frequency range. Within this band, one of the most promising sources are double white dwarf (DWD) bi- 
nary systems, as it is expect ed that a large fraction of main-sequence binaries end their lives as close DWDs 



(Iben fc Tutukov 



1984, 



19861 ). For this reason, the GWs emitted by these systems in our Galaxy may form a 
background noise in the low frequency (< 3 x 10^'^ Hz) band of LISA. The population of DWDs in our Galaxy 
is expected to be dominated by systems that undergo two distinct, long-lived phases of evolution: an "inspiral" 
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phase, where both the stars are detached from their Roche lobes and the loss of angular momentum in the form 
of GW emission causes the two stars to slowly spiral in towards each other; and a "stable mass transfer" phase, 
where the less massive star fills its Roche lobe initially and starts transferring mass steadily to its c ompanion. 
An e xample of the stable mass transferring systems are the AM CVn type systems, of which 18 (jNelemans 
2005 ) are known through electromagnetic observation^. 



Apart from DWDs, neutron star white dwarf (NSWD) binary systems are also one of t he promising sources 
of GWs for LISA. Various authors (|Kim et al.ll2004 ICooravll2004t iNelemans et al.l 120011) have estimated the 
GW background from these systems and concluded that the number of NSWD systems detectable with LISA is 
1-2 orders of magnitude less than DWD systems. Similar to DWDs, NS WD systems also undergo inspiral and 
stable mass transfer phases. S pecifically, several studies (see for example: iNelson et al.lll986l : lBelczvnski fc Taam 



2004a 



Nelemans et al.ll2006l ) have suggested that one of the possible formation scenarios of the so-called ultra 



compact x-ray binary (UCXB) systems, with orbital periods < 80 minutes, is that a low mass white-dwarf 
donor (< O.IMq) transferring mass to an accreting neutron star (NS) primary in a short orbit. In this scenario, 
a detached NSWD binary system initially evolves to a minimum orbital period as angular momentum is lost 
from the system due to GW radiation. At this minimum orbital period, the companion white dwarf (WD) star 
starts filling its Roche lobe and transfers mass to the NS and the system evolves to longer orbital periods. At 
present, there are 12 known UGXB systems with measured orbital periods, see Tabl({TJ Som e of these systems 
have a ccreting millisecond pulsars (XTE J1807-294, XTE J1751-305 & XTE J0929-314: see iMarkwardt et al 



(|2002f ) , for example ) and several of them are found in globular clusters as the stellar density and close encounters 
are more common (lGlarklll97 5l:E va nova et al.l 120071 ). The total pop ulation of field UCXBs may be low, ~ 10 



( Belczvnski &: Taam 2004a[" Coorav 2004) and theoretical studies by Belczvnski fc Taam ( 2004b ) indicates that 
even at the Galactic center, accreting NS systems do not contribute much to the faint x-ray population. 

In general, t he capabilities of LISA as a G W detector are usually discussed in the context of the log (h)- 
log (/) domain. iKopparapu fc Tohlind (j2007l ) proposed that an analogy can be drawn between the astron- 
omy community's familiar color-magnitude (CM) diagram and LISA'S amplitude-frequency diagram. For LISA 
sources, an analogous quantity to absolute magnitude M is log(r/i), where r is the distance to the source. The 
underlying physical properties of compact binary systems such as DWDs and NSWDs, their evolution, and 
their relationship to one another in the context of stellar populations can be ascertained only if the observa- 
tional propert ies of such systems are displa yed in a \og{rh) — log(/) diagram, rather than in a plot of log(/i) 
versus log(/). Kopparapu fc Tohlind (2007) discussed the DWD binary systems in this context of "absolute" 
amplitude-frequency domain assuming that the donors in these systems follow zero temperature mass-radius re- 



l ation . Here , we will consider the effect of "warm" donors as proposed in lDelove fc BildstenI (j2003t ): lDelove et al 



(|2005l . 120071 ) and extend the discussion to NSWD binary systems and compare the population boundaries be- 
tween both of them. We also discuss the limits and applicability of these population boundaries in the log(rft,) 
- log(/) space within the context of observed AM CVn and UCXB systems. 



2. Evolution of WD binaries in the amplitude-frequency domain 

As mentioned in the introduction, since the donor is a WD star and starts filling its Roche lobe to begin mass 
transfer phase, we will denote the donor WD with subscript d and the accreting companion as a. This notation 
will be followed even during the detached inspiral phase of the evolution. Hence, the total mass Mtot = J^Id + M^. 
and the mass ratio q ~ AI^/Mg_. Also we assume that the maximum mass of a WD to be Chandrasekhar mass. 



^Two controversial candidate sy s tems, RX J0806+15 and V407 Vul may change their number to 16. See lCropper et al 
IWu et al.l l l2002h : lMarsh fc Steeghj 1 I2OO2I) for more details. 
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Mch = IAAMq, the minimum mass of a NS is 1.2il/0 and the maximum mass of NS to be SAfg 



In the case of a detached DWD or a NSWD system inspiraUing as a resul t of lo ss of angular momentum 
due to gravitational radiation, we can write, as shown in lKopparapu fc Tohlind (j2007l ). as 



rh„ 
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where the dimensionless mass parameter. 
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M. = MtotQ'^^^ is the chirp mass and Q = q/{l + q)^. Notice that, for DWDs, the maximum value of K = 1 
occurs at M = 1.25Mq {q = 1, Mtot = 2.88Afo) but for NSWD binaries, the maximum value of K is 1.39 and 
occurs at = 1.75Af0 {q = 0.5, A'/tot = 4.32A'/q). This is purely due to the upper limit on the mass of the WD, 
Afch. Similar to a DWD system, where a _fir = 1 creates a maximum inspiral trajectory in the log[r/inorm]— log[/] 
plane, a, K = 1.39 makes the limiting inspiral trajectory for NSWD binary systems in this plane. This limiting 
trajectory, for both DWD & NSWD is plotted in Figs. [T]&[2l respectively, as a red line with a slope of 2/3 (see 
Eq.©) beyond which a detached inspiralling DWD or NSWD system can not be found. 

The phase of a detached inspiral evolution terminates when the low mass WD companion comes into contact 
with it's Roche lobe to initiate a phase of mass-transfer. Since the donor is a WD star this contact period can 
be found by equating the radius of the star with the Roche lobe radius. He re, for the radius of the W D, we 
use the mass-radius (M-R) relationship of non-zero entropy donor models of Delove fc Bildsten ( 2003h. These 



models were initially developed for ultra-compact x-ray binaries (UCXB) with WD donors. Later iDelove et al 



([20051) applied them to model the donors in AM CVn systems formed through WD channel. For these models 
the donors are assumed to be fully convective and corresponding donor evolution is consi dered to be ad i abatic 
in nature due to the large mass- transfer rates produced in AMCVn systems. However, iDelove et al.l (j2007() 



showed that the assumption of adiabatic evolution is applicable during the beginning of mass-transfer phases 
(when the mass-transfer rates are high) and may not work when the system evolves to larger orbital periods 
with low mass-transfer rates. Instead, they have identified new evolutionary phases and we will also discuss the 
impact of this most recent study on our population boundaries. We illustrate these boundaries by assuming 
that the donor composition is He, but it is straight forward to extend the discussion to carbon (C) & oxygen 
(O) donors. In this section, we will discuss in detail the population boundaries of both DWD & NSWD. In i|3] 
we will discuss the applicability of th ese boundarie s with respect to LISA observations, along with the effects 
of new phases of evolution detailed in IDelove et al. ( 2007 ). Finally, we summarize in [jH 



2.1. Population boundaries for DWD systems 

In Fig. [Tl We show the population boundaries for both inspiralling & mass-transferring phases of DWD 
systems. The (red) line with slope = 2/3 shows the maximum inspiral trajectory {q = 1, A^tot = 2.88Af0) from 
a detached DWD system. Above this boundary, no DWD binary can exist. The lower (green cross) curve shows 
the locus of the termination (contact) points of all detached inspiralling DWD binary systems, which have q = 1, 
assuming that the WD is a He donor and has a temperature of T = W^K (which is equivalent to the boundary 
obtained with T = M-R relation; along this isotherm, Aftot decreases from top to bottom). For comparison, 
we have also shown the contact boundary for q = I, ge nerated by Eggleton's mass-r adius (M-R) relationship 
for T = white dwarfs, (blue dashed line), as quoted by Verbunt fc Rappaport ( 19881 ) and also by Marsh et al 



(j2004[ ). At higher amplitudes and frequencies, both the curves overlap each other. This is due to the overlap 
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of M-R relations at higher masses as these objects are supported by degeneracy pressure and the T = M-R 
relation that we are using is applicable for fully degenerate He donors. At lower amplitude and frequencies, the 
contact boundaries diverge slightly. This is because there is a turn over in the M-R relation at low masses for 
T = W^K as the Coulomb interactions dominate thermal contribution at that low temperatures and the M-R 
relations deviate. So this (g = 1) contact curve represents the boundary beyond which no DWDs with He donors 
will be found. The corresponding boundaries for C & O would lie just below this contact boundary, except that 
all three (He, C & O) boundaries merge at high frequencies and deviate at low frequencies. 

Once the mass transfer phase is initiated, the orbital separation a starts increasing and the system evolves 
to lower amplitudes and frequencies. The evolutionary trajectories can be then traced out by assuming that 
the donor WD star is marginally in contact with it's Roche lobe and conservative mass transfer (CMT) holds 
true. For DWDs, Since we assume that the maximum mass in a DWD system is Mch, we can generate a upper 
boundary beyond which no mass-transferring DWD system can exist. This limiting boundary for DWDs can 
be obtained for systems with Aftot > Mch (because, then, eventually Ma will exceed Mch) when their q drops 
below a critical ^ch 

Mtot , . 

Accordingly, the (green) dashed line in Fig. [T] indicates an isotherm with Log[T] = 7.5 and = 1.44M0. We 
note that Log[r] = 7.5 upper boundary is likely a generous overestimate of the available phase space even if the 
donors were to evolve adiabatically from contact. Assuming that there are no donor WDs in a DWD system 
with Log fr] > 7.5, this curve rep resents the boundary beyond which no mass transferring DWDs can exist. 



Since the lDelove fc BildstenI (j2003[ ) donor models have two branch nature, with a minimum attainable mass, it 
is reflected in this boundary as it curves up, with a minimum attainable GW amplitude for a given isotherm. 
Corresponding boundary for donors with C & O composition would lie below this He curve because a higher 
temperature is required for a C/0 donor to fill the Roche lobe than for a He donor. For comparison, a similar 
boundary that arises from T = M-R relation is also shown as (light blue) dot-dashed curve. As is the case 
with q = I contact boundary, these two curves overlap at high frequencies. 

From the above discussion, we note that all mass-transferring DWDs with different donor masses, com- 
position an d temperatures should be co nstrained within these two boundarie^. Compared to T=0 model, as 



assumed in lKopparapu fc Tohlind (|2007l ). hot donors in mass-transferring DWDs occupy larger fraction of space 
in rh~ f domain. This is because hot donors are also more massive, therefore have higher intrinsic GW amplitude 
and come into contact at lower frequencies because of larger radius. 

The inspiral phase of a WD binary evolution is driven by the loss of the angular momentum due to GW 
radiation and the inspiral evolutionary time scale Tchirp can be written as: 



q 



647r2 \ rhnoimP 



(4) 



It is worth noting that the t ime scale of evolution in both inspiral and CMT phases is of the order of ~ Tch 



irp 



(jKopparapu fc Tohlinell2007l ). Therefore we have drawn "chirp" isochrones also in FigUand each isochrone has 
a slope of —2 in this log(r/inorm) - log(/) space, as can be seen from Eq.Q due to the dependence on the product 
of r/inorm and This implies that a given WD binary system, in either inspiral or CMT phase, spends Tchirp 
amount of time corresponding to the GW frequency that the binary is emitting. At lower frequencies, the value 
of Tchirp is larger compared to that of at higher frequencies. A consequence of this behavior is that more binary 
systems accumulate at lower frequencies. As mentioned in the introduction, the background noise arising due 



For C & O composition donors, as discussed above, there is a slight deviation at the lower contact boundary. 
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to millions of DWDs at the lower frequency band (< 3 x 10 Hz) of LISA is due to the fact that the value of 
Tchirp for these systems is 10^° years, approaching Hubble time (see Figl^]). 




"q=?1 inspiral 

q=1, Log[T]=4 contact 

Log[T]=7.5 contact 

T=0 contact. 

Ma = M(,h,T = 0''-, 

AMCVn 
I 



-2.5 -2 -1.5 

Log[f(Hz)] 

Fig. 1. — Population boundaries for DWD systems. The (red) line with positive slope represents the maximum 
inspiral evolution trajectory for DWD systems, beyond which none will be found. This happens at q = 1, ilftot = 
2.88Mf;^. The bottom (gre en) curve with crosses is drawn using the M-R relation of isentropic donor models from 
Delove fc BildstenI |2003h . It represents the termination (contact) boundary for q = 1 systems, assuming that 
the donors are of He composition with T = W^K. For comparison, contact boundary derived using Eggleton's 
T = M-R relation for fully degenerate He donors is also shown as (blue) dashed line. The top (green) dashed 
curve represents the boundary at which 7\/a = Mch, assuming that donors in DWD systems would have a 
maximum T = No mass-transferring DWD would be found above this boundary. For comparison, the 

(blue) dot-dashed curve shows a similar boundary for T = donor models. Also shown are constant Tchirp lines 
along with observed AMCVn systems (blue dots; see text). Explanation for vertical (brown) dot-dashed line is 
given in ^ 
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2.2. Population boundaries for NSWD systems 

In Fig. [21 population boundaries for both NSWD systems is shown in the log(r/inorm) - log(/) space. The top 
(red) curve represents the maximum inspiral boundary for detached NSWD binary systems (K = 1.39) beyond 
whic h none will be found. The bottom (green) curve with crosses represents the lower contact boundary, drawn 



using lDelove fc BildstenI ([20031) hot donor models, assuming that the donor is of He composition with T — W^K. 
The mass of the NS on this boundary is set to the minimum mass assumed, 1.2Mq. For comparison, a similar 
boundary using T = M-R relation for He donors is also shown as (blue) dashed curve. As is the case with 
DWDs, these two boundaries overlap at higher amplitudes and frequencies because at higher donor masses, the 
two M-R relations match as they are supported by degeneracy pressure. Furthermore, there is an overlap region 
between the lower NSWD contact boundary (green dashed curve in Fig. ^ & the Ma = M^h DWD contact 
boundary (green dashed curve in Fig.[T]), indicating both types of systems can exist within this region. Similarly, 
an upper contact boundary can also be drawn for NSWDs, shown as (green) dashed curve in Fig. [21 assuming 
the maximum mass of the NS is 3.0Mq and that the maximum donor temperature in a NSWD system to be 
T = W-^K. The light blue dot-dashed curve is the upper contact boundary for T = donors. Just like DWDs, 
for NSWD systems with Mtot > 5Mq, the critical mass ratio (/ns below which a NSWD binary can not exist, is 

'^^^ = 3Mo-^ 

So all the mass-transferring NSWD binary systems are bounded by the upper & lower (green) curved and all 
the detached inspiralling NSWDs will be bounded between the top (red) curve and the bottom (green) curve. 
Fig. [2] also shows chirp isochrones, indicating the evolutionary timescalcs. 



*Thc lower boundary is applicable for He donors. 




-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 



Log[f(Hz)] 

Fig. 2. — Population boundaries for NSWDs in the log[r/inorm]-log[/] plane. The red line represents the maximum 
inspiral trajectory from a detached NSWD binary system {q = 0.5, Afmtot = 4.32M0). The bottom curve with 
(green) crosses represents the contact boundary, assuming minimum mass of the NS to be 1.2Mq with He donors 
of T = 10^ A'. The blue dashed line is drawn assuming a T = M-R relation for fully degenerate He donors. 
These two boundaries overlap at high masses (high amplitude & frequency) and deviate slightly at lower masses. 
The top (green) dashed curve shows the upper bound for mass-transferring NSWD systems, assuming that the 
maximum donor temperature in a NSWD system to be T = lO'^-^K and maximum mass of NS as 3.0 Mq. To 
compare, an upper boundary drawn using T = donor model is also shown as (light blue) dot-dashed curve. 
Also shown are currently known UCXB systems (blue dots) along with chirp isochrones. 
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Discussion 



The population boundarie s for DWD & NSWD binary systems discussed in previous sections were generated 



using iDelove fc Bildstenl (j2003f ) non-zero entropy donor models, which assume that the donors are fully convec- 
tive and undergo adiabatic evolution throughout the mass- loss phase. Recently, iDelove et alj (|2007l ) showed that 
these assumptions may not properly estimate the donor's orbital period evolution and specifically the donor's 
adiabatic evolution may not hold true for whole mass-transfer phase. Instead they identified three distinct 
evolutionary phases and here we discuss the implications of their new findings on our population boundaries. 

During the first phase, which happens during the mass-transfer "turn-on" phase (when the donor comes 
into contact initially), the radius of the donor decreases , the mass-transfer rate Md increases and the orbital 
period Porb continues to de crease until the don or radius reaches it's minimum value (Md becomes maximum) 
and starts expanding again. IDelove et al.l (|2007l ) calculated that this turn on phase lasts up to ~ 10^ years. In 
the second phase, the donor responds to the mass loss adiabatically and starts expanding, which is considered 
to be the normal AM CVn phase. But this phase of adiabatic expansion ends and a third phase of evolution 
begins at around Poib ~ 45 min, when the mass-transfer rate (and the donors thermal time) drop enough for 
the donor to cool and start contracting to a fully degenerate configuration, stalling the Porb evolution. 

In the case of DWDs, the initial turn-on phase will result in the q = I conta ct boundary (green curve wit h 
crosses in Fig. [ij to have a width, instead of a sharp boundary as discussed in iKopparapu fc Tohlinel |2007l ). 
This is because it is calculated assuming that once the system comes into contact, it would evolve towards lower 
amplitudes and frequencies, whereas Porb decreases in the turn on phase, to a minimum even after the initial 
contact. Accordingly, there will be a slight overlap of lower contact boundaries for He, C & O donors. Once the 
system evolves off this initial contact phase, the donor expands adiabatically in response to the mass-transfer 
and the system follows a typical AM CVn evolution, where the GW amplitude & frequency keeps decreasing 
as Porb increases. Assuming an adiabatic evolution means that the cooling time of the donor is longer than 
the mass-transfer time-scale (Md/Md), it will in turn affects the orbital evolution time-scale. Accordingly, the 
do nor follows a trajectory wher e it passes through different isotherms of decreasing q, after the initial contact. 
As lKopparapu fc Tohlinel ([20071) illustrate, the lower right curve in Fig. [T] shows the contact boundary for g = 1 
systems and similar contact boundaries for lower g's would lie to the left of it, but the contact boundaries will 
shift towards lower amplitudes and frequencies. This phase of adiabatic evolution comes to an end between 
Porb ~ 40 — 55 min, when the donor starts to cool and contract e ventu ally towards a fully degenerate star. This 
will drastically (almost an order of magnitude, see IDelove et al.l (j2007l ) Fig. 15) reduce Md at these long orbital 
periods. Hence, after this range of Porb, the systems GW frequency evolution slows down in accordance with 
the drop in Porb evolution. If the donor has cooled enough to approximate it as a P = degenerate model, 
then it may lie close to one of the T = contact boundaries, depending upon its q. But these P = boundaries 
are bounded within the region constrained by 7\/a = Mch (light blue dot-dashed) curve and the q = 1 contact 
boundary (green curve with crosses in Fig. [T]), because these boundaries are drawn assuming that the donor is a 
fully degenerate star with P = and such a system can not exist beyond these curves. Therefore, we have drawn 
a vertical (brown) dot-dashed line in Fig.[T]at Porb = 55 min (Log/ = 3.21) beyond which we would expect these 
systems to be within the P = contact boundaries. Note that this vertical line is not a sharp boundary: some 
systems may not reach a fully degenerate configuration by this Porb- Rather, a system's Porb evolution slows 
down at a particular GW frequency once they start cooling towards a degenerate configuration at the above 
ment ioned Porb. Fig. [T] al so s hows the obse r ved A M CVn type systems, for which the masses and Porb are taken 
from belove et aP (2005) & Roelofs et al. ( 2007 ). Couple of them are clearly outside P = region; for some 



of them, the system's mass function limits the minimum donor mass to a value above that of a Roche-filling 
P = donor at the same Porb- It is very difficult to know their exact temperature and/or composition purely 
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from GW observations, as there is a good overlap of systems with these characteristicfU. 

For NSWDs, similar to DWDs. the contact boundaries with He composition drawn assuming NS mass is 
I.2M0 (green curve with crosses ) and 3.0Mq (green dashed curve) will also have a width due to decrease in 
Porb even after the contact. Furthermore, the lower contact boundary lies below the upper contact boundary 
for DWDs (Ma = IAAMq). This will make the contact boundaries to "overflow" into the DWD region and 
LISA observations may not be able to distinguish these two types of systems in this region. Here also the 
donor undergoes adiabatic evolution after reaches a maximum and consequently enters the cooling phase at 
Poib ~ 55 min. Accordingly, in Fig. [21 the brown line shows the orbital period of this transition phase. Beyond 
this line, the donors should start cooling and GW frequency slows down accordingly. Fig. [2] also shows the 
currently observed UCXBs, which we assume to be mass-transferring NSWD binary systems. Table [1] gives the 
values of the masses of donors and orbital periods that we used. Some of the UCXBs shown in Fig. [2] lie below 
the lower contact boundary (green cross curve), indicating that probably the donors are not of He composition. 
Although, this contact boundary has a width, and hence they may have He composition, it is unlikely that 
we can observe them during this relatively short lived phase. Moreover, these systems are plotted assuming 
the minimum mass of the donors mass range derived from observations, so this provides additional uncertainty 
in determining the composition of donors. It is unlikely that LISA will be able to observe cooling donors in 
either DWD or NSWD systems because of the instrumental and/or DWD background noise. In FiglHl we plot 
the known UCXBs and AM CVns on top of LISA's sensitivity curvcj^ (SNR = 1) to assess the detectability of 
these systems. In the case of known UCXBs, it is clear that only one system (4U 1820-30) has enough signal 
strength to be visible to LISA, whereas some of the known AM CVns emit GWs above the instrumental and 
DWD background noise. 

Recalling from IJI] in order to transform from log(/iiiorm) to log(r/inorm) space, we need to know the distance 
r to the binary system. The relati on between the unknown bina ry parameters r, Mtot and q and the observables 
hnoim, f and / can be written as ( Kopparapu fc Tohlinel 2007 ) 



ML 



267r2G5 p 



"norm 



(1 + 9) 

j^^ii ii'normj 

where g ~ in the inspiral phase of evolution and hence, it is easy to determine r from /inorm, / and / through 
Eq.(I7]). For mass-transferring systems g is a function of Mtot and q, and they can be related to the observable 
/ by the requirement that in the mass transfer phase, i?d = i?L- 

The determination of r and/or the masses of the stars in DWD/NSWD binary system depends on the 
determination of /. If an / can not be measured for a system, then there is no way to tell whether that system 
is a NSWD or DWD system based only on LISA observations. If an / can be measured, and if it turns out to 
be negative, then it is possible that particular system is a mass-transferring system (DWD or NSWD). But as 
shown in Figs.[l]&[2l it will still not be possible to know the type of the system, at least for low frequency sources 



^Out of these systems, CE 315 (the system with lowest GW amplitude and frequency in Fig. [TJ has Po^b = 65 min and lies to 
the left of the Porb = 55 line and below T = 0, Ma = boundary. This may seem to im ply that this system may be cooling off 
and trying to reac h a T = degenerate configuration. But in Fig 1. of lDelove et~al] (|2005|) . they give a tem perature range of Log 
T = 4.0 — 6.5 and iBildsten et al. I 1 I2OO6I) noted that this system may have a hot donor. If that is the case, iBildsten et aU l l2006h 
note that this donor may have evolved with constant entropy that it was b orn with or may have been heated by either the disk or 
the accretor. However, if irradiation is the cause, then belove et al. note that it will delay the onset of donor's cooling and 

increases the temperature of the donor. So though this system lies below the T = 0, A/a = Afch boundary, it does not necessarily 
mean that the donor can be approximated as a zero-temperature object. 

®http:/ /www. srl. caltcch.edu/ shanc/scnsitivity/MakcCurvc.html 
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(Log[/] ^ — 2). There is a fairly good overlap in / between DWD & NSWD systems in this region because of 
the non-zero entropy nature of the donors and also due to the lower limit on the mass of the NS (1.2Afo). But 
for high frequency mass-transferring sources (Log[/] ^ ~2), it may still be possible to know the type of the 
system, as the overlap region reducef0. The same thing can be said about inspiralling systems because there is 
a large overlap of DWD and NSWD inspirals at low frequencies and even a measurement of positive / would 
unlikely be able distinguish these two types of inspiralling systems. 



4. Summary 

In the previous sections, we have outlined the construction of population boundaries to illustrate the 
various evolutionary phases that a DWD or a NSWD binary system would undergo in th e distance independent 



"abso lute" amplitude-frequency domain (log[r/inorm] - log[/]) of LISA. In an update to iKopparapu fc Tohline 



(|2007l ). who assumed the donors to be fully degenerate T = He stars, we consider that the dono rs in these 



systems follow the M-R relationship of non-zero entropy donor models of iDelove fc BildstenI ([20031) assuming 
He composition. These models assume fully convective and adiabatically evolving donors during the whole 
episode of the mass-transfer phase. Figs. [T] & [2] show that these "hot" donors occupy a larger fraction of the 
rh — f space, than the T ~ donors, because hot donors are also more massive, increasing their intrinsic GW 
amplitude. At high frequencies, both the models match each other because at these high masses, the hot donors 
are supported by degener acy pressure and the M-R relations match. At low frequencies, the T = model and 
Delove fc BildstenI (j2003l ) T = lO^if full model (which is equivalent to T = 0) diverge slightly because the 



Coulomb interactions dominate thermal contribution and M-R relations deviate. 



We also discussed the implications of new evolutionary phases found by iDelove et al.l (|2007l ) on our popu- 
lation boundaries. The initial "turn on" phase, where reaches it's maximum value and Poib decreases even 
after contact, will result in contact boundaries having a width, instead of an abrupt cut-off. This will cause 
some overlap onto the C fc O isotherms, which lie below He boundary. Soon after the system evolves from this 
initial contact, the second phase starts where the donor undergoes adiabatic evolution in response to the mass- 
transfer and the system will follow a typical AM CVn evolutionary trajectory with decreasing GW amplitude 
and frequency. This will continue until the donor begins to cool and the radius starts contracting. Accordingly, 
it will reduce Md, and Porb (and GW frequency) evolution slows down and may stall once the donor reaches 
fully degenerate configuration. But since LISA's sensitivity in this range is limited by instrumental and DWD 
background noise, this part of the evolution (or cooling donors) will probably be not observable by LISA. 

It is unlikely that determination of / will shed light on the type (DWD/NSWD) of the low frequency 
systems without the help of independent electromagnetic observations. This is because there is fairly a good 
overlap of NSWD and DWD systems within the resolvable frequency regime (Log[/] > 3) of LISA Q But for 
high frequency sources, it might be possible to distinguish between them. Combined with the expectation that 
the relative population of NSWD is low compared to DWDs and short period systems do not stay longer at 
those periods, it is likely that LISA will be able to measure / for more number of DWDs than NSWDs. It 
will be interesting to see how the high frequency regions constrained by these boundaries would be populated 
through LISA observations and whether indeed we will know the nature of these systems. 



''There still will be some uncertainty for systems with NS mass lower than M^y^, but for NS masses higher than A/^jj, LISA 
should be able to distinguish both types of systems through the measurement of /. 

® Since this region is occupied by both inspiralling & and mass-transferring systems, it may happen that more than one source 
can reside in a frequency bin, but a measurement of / would probably indicate the evolutionary phase. 
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Table 1. Observed and derived properties of some of the know n UCXBs (|in 't Zand et al j[ 2007f) shown in 
FigH Md indicates minimum mass of the d onor. References: (l)lTarana et all (l2007h (lA l lCumminj |2003h 
(2) Iwang fc Chakrabartvl tooi) (2A) IWang fc Chakrabartvl (l2Q04h O l lSidoh et all I2OO6I) (3A) 



Papittol (l2008h rsllCallowav et all tooi) (8A ) 



Gallowav et al 



( 200: 



I(10A) 



Chakrabartyr(| 199 



2OO5I) f4A) i McNamara et alJ (l2004h f5) IPalanga et all j2005l) 
J2OO7I) ('BATlChakrabartv l dlQQSh ('Tlcierhnskij tO()^ (7A) 

i) (llA) 



Nelemans et alJ (|2004l ) (12) iKrimm et alJ (|2007l ) (12A) http: / /web.mit .edu /newsoffice/2007/pulsar-09 IS.html 



Name 


Orbital period 
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M@ 


Reference 
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Reference 




(min) 


for mass 


(kpc) 
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Fig. 3. — The known UCXBs (blue stars) and AM CVn systems (light blue stars) plotted on top of LISA's 
sensitivity curve for SNR = 1. In case of UCXBs, only one system (4U 1820-30) has enough signal strength to be 
visible to LISA, whereas couple of known AM CVns emit GWs above the instrumental and DWD background 
noise (green curve). 
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